Confronting quasi-exponential inflation with WMAP seven 



(N 
O 



o 
u 

Oh! 

o 



> 
m 

o 



Barun Kumar Pal^*, Supratik Pal^^^^' and B. Basu^^ 
Physics and Applied Mathematics Unit, Indian Statistical Institute, 203 B.T. Road, Kolkata 700 1 08, India 
^ Bethe Center for Theoretical Physics and Physikalisches Institut der Universitdt Bonn, Nussallee 12, 53115 Bonn, Germany 

We confront quasi-exponential models of inflation with WMAP seven years dataset using Hamilton 
Jacobi formalism. With a phenomenological Hubble parameter, representing quasi exponential 
inflation, we develop the formalism and subject the analysis to confrontation with WMAP seven 
using the publicly available code CAMB. The observable parameters are found to fair extremely well 
with WMAP seven. We also obtain a ratio of tensor to scalar amplitudes which may be detectable 
in PLANCK. 



I. INTRODUCTION 

Present-day cosmology is inflowing into an era where 
it is becoming more and more possible to constrain the 
models of the early universe by precise data coming from 
highly sophisticated observational probes like WMAP [l| , 
PLANCK 0, SDSS [|, ACBAR [J, QUaD 1]. Such 
observations are gradually leading theoretical cosmology 
towards the details of physics at very high energies, and 
the possibility of testing some of the speculative ideas 
of recent years. Inflation - the most fascinating among 
them - was first proposed back in 1981 by Guth 0, and 
has since been developed by a wide variety of probable 
models (rl-fl^. However, the inflationary scenario is more 
like a paradigm than a theory, since a specific compelling 
model is yet to be separated out from the spectrum of 
possible models and alternatives. 

So far the most appealing prediction of inflationary 
paradigm was that it may generate spectra of both den- 
sity perturbations and gravitational waves. Before the 
detection of CMB anisotropies by COBE [l^ , cosmolog- 
ical observations had limited range of scales to access, 
and it was sufficient to predict from an inflationary sce- 
nario a scale-invariant spectrum of density perturbations 
and a negligible amphtude of gravitational waves. How- 
ever, since COBE ^ and, of late, WMAP Q came for- 
ward, the spectra are now well-constrained over a wide 
range of scales ranging from 1 Mpc upto 10,000 Mpc. So 
the inflationary predictions should now be very precise to 
incorporate latest observations, say, from WMAP seven 
year run. The ongoing satellite mission PLANCK 01 will 
restrict the theoretical predictions further. It promises to 
survey the ratio of the tensor to scalar amplitudes r upto 
the order of 10"^ [16] (orders below WMAP [1] predic- 
tions r < 0.36 at 95% Confidence Level [C.L.] [17]) so as 
to comment more precisely on primordial gravitational 
waves and the energy scale of inflation. PLANCK may 
even discriminate single field infiationary models from 
multi-field models by detecting (provided /atl > 5) large 
non gaussianity. So this is high time to confront dif- 



ferent class of inflationary models with latest data and 
forthcoming predictions. 

The usual technique used in the inflationary scenario 
to solve the dynamical equations is the slow roll approx- 
imations [18|. But it is not the only possibility for suc- 
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cessfully implementing models of inflation and solutions 
outside the slow roll approximations have been found [13] ■ 
To incorporate all the models irrespective of slow roll ap- 
proximations Hamilton Jacobi formalism [21rf23i | turned 
out to be very useful. The formulation is imitative by 
considering the inflaton field itself to be the evolution 
parameter. The key advantage of this formalism is that 
here we only need the Hubble parameter H{(f) to be spec- 
ified rather than the inflaton potential V{(j)). Since H is 
a geometric quantity, unlike V , inflation is more natu- 
rally described in this language 0, Further, being 
first order in nature, these equations are easily tractable 
to explore the underlying physics. Once H{<j)) has been 
specified, one can, in principle, derive a relation between 
(f) and t which will enable him to get hold of H{t) and 
the scale factor a(t) therefrom. Therefore Hamilton Ja- 
cobi formalism provides us with a straightforward way of 
exploring infiationary scenario and related observational 
aspects. 

In this article, we intend to confront quasi-exponential 
infiationary models with WMAP seven years data [l| 
with a phenomenological Hubble parameter following 
Hamilton Jacobi formalism. The absence of time depen- 
dence in the Hubble parameter produces de-Sitter infia- 
tion which was very appealing from theoretical point of 
view but its acceptability is more or less limited con- 
sidering present day observations. So, certain deviation 
from an exact exponential inflation turns out to be a 
good move so as to go along with latest as well as forth- 
coming data. In general these models are called quasi- 
exponential infiation. Our primary intention here is to 
confront this quasi-exponential inflation with WMAP 
seven using the publicly available code CAMB 24| . Nev- 
ertheless, as it will turn out, the analysis also predicts a 
detectable tensor to scalar ratio so as to reflect signifi- 
cant features of quasi-exponential inflation in forthcom- 
ing PLANCK data sets as weU. Thus, PLANCK may 
directly verify (or put further constraint to) our analysis 
by detecting (or pushing further below the upper bound 
of) gravitational waves. 
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II. MODELING QUASI-EXPONENTIAL 
INFLATION BY HAMILTON JACOBI 

The Hamilton Jacobi formalism allows us to express 
the Friedmann equation in terms of a first order differ- 
ential equation ^ 



H {ct>) 



2Ml 



1 



1/(0) 



and the evolution of the scalar field 
0= -2MlH'{(t)) 



(2.1) 



(2.2) 



where a prime denotes a derivative with respect to the 
scalar field </> and a dot a time derivative. The above two 
equations govern the infiationary dynamics in Hamilton 
Jacobi formalism. The shape of the associated potential 
can be obtained by rearranging terms of Eqn. (|2.1[) to give 



F(0) = 3M|,ij2(0) 
where tu has been defined as 



1 

1 - -€h 

3 " 



tH = 2Mf, — 




(2.3) 



(2.4) 



We further have 



: H{<l>f [1 - en] 



(2.5) 

So < 1 implies accelerated expansion. It is customary 
to define another parameter by 



riH = 2M^ 



H 
IT 



(2.6) 



It is worthwhile to mention here that the parameters en 
and rjH are not the usual slow roll parameters. But in 
the slow roll limit ejj — )■ e and i]h ^ V ~ ^ 1 ^ ^^'^ V 
being usual slow roll parameters. 

Let us consider a phenomenological Hubble parameter 
representing quasi exponential infiation 



i7(0) = H„if exp 



Mp 



P{1 + 



- 
Mp ' 



(2.7) 



where p is a dimensionless parameter and is a con- 
stant having dimension of Planck Mass. The value of the 
constants can be fixed from the conditions for successful 
inflation and the observational bounds. As it will turn 
out in subsequent analysis this Hubble parameter can in- 
deed be cast into a form of quasi exponential infiation for 
some choice in the parameter space. 

The two parameters tn and r\H in the Hamilton Jacobi 
formalism take the form 



2-1 



2P + 2P^ 



Mp . 



p2(l. 



Mp 



(2.8) 
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FIG. 1: Variation of the fourth root of the potential with the 
scalar field 
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FIG. 2: Variation of the equation of state parameter with the 
scalar field </!>, the horizontal line represents uj = 



Now, the condition for inflation may be put forward 
using equation of state parameter which has the form 



_ P{4>) 



1 



(2.9) 



P(0) 3p2(l + Ji-Y 
From the constraint w < — ^ during inflation, we get 



Mp ■ 



ever, at the end of inflation a; > — ^ allows us to estimate 



So for almost any value of p infiation is achievable. How- 
ever, at the end of inflatic 
an upper bound for p as 

V2 



\V\ 



< 



(2.11) 



where ipend is the value of the inflaton field at the end 
of inflation. In order to implement a successful model of 
inflation both the restrictions on lj should be satisfled. 
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7.448 


70 
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TABLE I; Different parameters for different values of p within 
its allowed range 



First note that if p is negative then from Eqn. (|2.2|) and 
Eqn. (l2.7p we have (p > 0, i.e. (j>{t) increases with time 
which incorporates the so called graceful exit problem (26l , 
[27| . We discard the values of p greater or equal to 1/2 
as well, since p > V2 would imply e// < 1 for any value 
of (f> giving rise to the same problem. As a result the 
feasible range for p would be: < p < \f2. Further, for 
sufficient inflation we also need \r\ii\ < 1 during inflation. 
And violation of that condition after inflation drags the 
inflaton towards its potential minima quickly. But if p < 
0.586 then \r\H\ will always be less than one. Imposing 
this condition further restricts the range of the otherwise 
free parameter p within 0.586 < p < -\/2. We shall use a 
representative value for the parameter p within the above 
range later on in this paper while confronting CAMB [23] 
outputs with WMAP seven 

Now to derive the expression for the scale factor we 
shall make use of the relation 



(2.12) 



Plugging Eqn. dsn]) into Eqn. ([22]) and using Eqn. ([2l^ 
we have 



a((/)) = fle exp 



(2.13) 



where i 



id) exp 



L)3 



a(0end) being the 



scale factor at the end of inflation. 

The primary quantities related to inflation have been 
summarized in Table ID We see that \t]h\ « 1 after the 
end of inflation, and so slow roll would be a very good ap- 
proximation throughout the inflationary period, though 
we do not utterly require them in this formalism [25j . 

That the above scale factor indeed represents quasi- 
exponential inflation will transpire from the following 
analysis. We first expand i?((/)) into the power series 
of 6 



H{4>) = Hi 



1 - 



1 



1 



1 



p(l + ^) 2!p2(i + _|_)2 



Mp ■ 



where we have defined Hi = exp 



(2.14) 
The above 



expansion, at next to leading order, gives rise to 



Mf 



P 



1 



(2.15) 



Here te = ^tend + (1 + Tif)^ ^^'^ *end corrcsponds 
to the time at the end of inflation. 

The combination of Eqn. ()2.13p and Eqn. ()2.15|) gives 
the scale factor as a function of cosmic time t 

a{t) aeexp[Hi{t-te)] (2.16) 

So the expression for the conformal time turns out to be 

1 



77 : 



Hia{t) 



(2.17) 



Thus, our analysis indeed deals with quasi-exponential 
inflation. The higher order terms in the expansion (|2.14p 
will, in principle, measure further corrections to the scale 
factor but the analytical solutions are neither always ob- 
tainable nor utterly required. Rather, one can directly 
confront the observable parameters with WMAP seven 
in order to constrain quasi-exponential inflation, as done 
in the rest of the article. 



III. PERTURBATIONS AND OBSERVABLE 
PARAMETERS 

A. Scalar curvature perturbation 

The fc*'' Fourier mode of the comoving curvature per- 
turbation, obtained by solving Mukhanov-Sasaki equa- 
tion in the slow-roll limit and adopting standard Bunch- 
Davies ;28] vacuum as initial condition, is approximately 
given by 

pgir; [|ln(ij'iae|r?|)j expHfc??] / _^ 

(3.1) 

So the dimensionless spectrum for the comoving curva- 
ture perturbation is. 



Pn{k) 



1 + k^TJ^ 



6 

-ln(i7iae|77| 
P 



4/3 



(3.2) 



Now to evaluate the spectrum at horizon exit k — aH, 
we notice that 



1 + k^ri^ 



6Hip^{t, ~ t) 



(3.3) 



To arrive at Ean. p.3p we have used Ean. (|2.17p and 
Eqn. (12.151) respectively. This reflects the features of 
quasi-exponential behavior where the effect of the evo- 
lution of the scalar field has been directly taken into con- 
siderations, without using k = —r]~^ a priori. We shall, 
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of course, use the relation k = —rj~^ finally while esti- 
mating parameters. Loosely speaking, the argument is 
analogous to putting a specific value for a variable after 
integration (which is more accurate) rather than putting 
it before (which may be a good approximation at the 
best). Somewhat similar treatment, though from a dif- 
ferent physical perspective, has been employed in [29]. 

Therefore the power spectrum for TZk evaluated at 
horizon crossing is now given by 
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-1.343 


-4.334 


-3.403 


0.878 



QA{k) 



P 



4/3 



'""^'A (3.4) 



P 



where we have defined A{k) = In (iJiagfc-^) . Also at 
horizon crossing d\nk = Hidt, so the scalar spectral in- 
dex looks 



TABLE IL Observable quantities as obtained from the theory 
of fluctuations 



Here tensor spectral index nx — and its running 
n'rp = 0, the results that we would obtain by the direct 
substitution of fc = ~7]~^ into Eqn. (j3.7l) . 

We are now in a position of dealing with one of the 
most important observable parameters, namely, the ratio 
of tensor to scalar amplitudes. In the present context we 
have found 



nsik) = 1 - 



4/3 



A{k) 



1/3 



_6_ 

p2 



16 



6A{k) 
P 



4/3 



6A{k) 



(3.5) 



[6yl(fc)]4/3p2/3 



(3.11) 



and its running 



n'sik) 



(|)'^'A(fc)V3_|>J^^(fc)l/3 + p 



In the next section we shall comment on its numerical 
estimates. 

The consistency relation can be obtained by combining 
Eqn. p.9| ) and Eqn. p.lip . which in this case turns out to 
be 
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'6A{k)' 








P 


) 


T = — 8nT 



Here note that we would get ng = \ 



(3.6) 



^1/4 



(3.12) 



33(fe) ^'^^ 



— g^^^p the expressions for the scalar spectral index and 

its running respectively if we had put k — ~r\^^ directly 
into Egn. iP^ . 



B. Tensor perturbation 

The power spectrum for the tensor modes representing 
primordial gravitational waves, obtained in a similar way 
as before, is given by 



Prik) = 2Ph{k) = 



Hi 



(1 -I- k^rf) (3.7) 



Now using Ean. (l3.3p we evaluate Eqn. (|3.7p at horizon 
crossing 

1 



PTik)\k=aH = 



2Hl 

^2 A/2 



1 



[6p2A(fc)]l/3 



(3.8) 



Corresponding spectral index can be derived immediately 
from Eqn.dSH) as 



nrik) = - 



1 



1 



(3.9) 



3 (6p2)l/3^(fc)4/3 _ ^(fc) 

and the running of the tensor spectral index is given by 



4[6p2A(fc)]i/3_3 



(6p2)l/3^(fc)4/3 _ A{k) 



(3.10) 



The result is slightly different from the known relation 
r = —Sht- In the literature there are other techniques 
for obtaining a modified consistency relation. Such a 
modified consistency relation can be found in any anal- 
ysis where higher order terms in the expansion of slow 
roll parameters are taken into account [30|. The consis- 
tency relation is also modified in the context of brane 
inflation [29l . [3ll | and non-standard models of inflation 
[3^ where generalized propagation speed (less than one) 
of the scalar field fluctuations relative to the homoge- 
neous background have been considered. Further devia- 
tion from the usual consistency relation can be found in 
(33} where tensor to scalar ratio has been shown to be 
a function of tensor spectral index, scalar spectral index 
and running of the tensor spectral index. Our approach 
is somewhat similar to these. 



IV. CONFRONTATION WITH WMAP SEVEN 

A. Direct numerical estimation 

In Table HI] we have estimated the observable parame- 
ters from the flrst principle of the theory of fluctuation as 
derived in the previous section, for two sets of values of p 
within its allowed range, 0.586 < p < \/2. For estimation 
we have taken = 2.27 x 10"^ Mp, = 7.5 x IQ-^^ 
and set the pivot scale at fcg = 0.002 Mpc^^. Table [H] 
reveals that the observable parameters as derived from 
our analysis are in excellent agreement with the current 
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TABLE III: Input parameters 
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Vo 
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13.708 


10.692 


0.2669 


0.7331 


0.0 


285.15 


14347.5 



TABLE IV: Different physical quantities as obtained from 
CAMB 



observations as given by WMAP seven years data for 
ACDM background [J. 

Nevertheless, present day observations are eagerly 
waiting to detect primordial gravitational waves. The 
current observational bound on the ratio of tensor to 
scalar amplitudes as given by WMAP is r < 0.36 at 
95% C.L. With the ongoing satellite experiment 

PLANCK promising to detect it down to the order of 
10^^ 16] in near future, this parameter is now playing a 
pivotal role (along with non-gaussianity) to discriminate 
among different models |34| . The numerical estimation 
reveals that the parameter r is of the order 10~^ for quasi- 



exponential inflation. Thus, PLANCK ^ may directly 
verify (or put further constraint to) our analysis by de- 
tecting (or pushing further below the upper bound of) 
gravitational waves. 



B. CAMB output and comparison 

In what follows we shall make use of the publicly avail- 
able code CAMB 24] in order to confront our results di- 
rectly with observational data. For CAMB, the Ean. (l3.4p 
has been set as initial power spectrum and the values of 
the initial parameters associated with inflation are taken 
from the Table HIl for p — 0.60 within its allowed range 
0.586 < p < V2. Also WMAP seven years dataset for 
ACDM background has been used in CAMB to obtain 
matter power spectrum and CMB angular power spec- 
trum. We have set the pivot scale at fco = 0.002 Mpc~^. 

Table Hn] shows input from the WMAP seven years 
dataset for ACDM background. 

Table IIVI shows the output obtained from CAMB, 
which is in fine concord with WMAP seven years data. 
The results obtained here are for a representative value 
of the parameter p — 0.6 within its allowed range. 

The curvature perturbation is generated due to the 
fluctuations in the inflaton and remains almost constant 
on the super Hubble scales. Long after the end of in- 
flation it makes horizon reentry creating matter density 
fluctuations through the gravitational attraction of the 
potential wells. These matter density fluctuations grew 
with time forming the structure in the universe. So the 
measurement of the matter power spectrum is very cru- 




FIG. 3: Variation of matter power spectrum PMat{k) with 
k/h m logarithmic scales for quasi-exponential inflation and 
the best fit spectra of WMAP7 for ACDM + TENS 



CMB TT Angular Power Spectrui 




FIG. 4: Variation of CMB angular power spectrum C; for 
quasi-exponential inflation, the best fit spectra of WMAP7 
for ACDM + TENS and WMAP7 data with the multipoles I 



cial as it is directly related to the formation of structure. 
In Figl3]the CAMB output for the variation of the spec- 
trum of the matter density fluctuations with the scale 
for quasi-exponential inflation and the best fit spectra of 
WMAP7 for ACDM + TENS [H have been shown and 
it represents true behavior indeed Q. 

In Figg] we confront CAMB output of CMB angu- 
lar power spectrum Cf" for quasi-exponential inflation 
with WMAP seven years data and the best flt spectra 
of WMAP7 for ACDM + TENS. On the large angular 
scales i.e. for low I, CMB anisotropy spectrum is dom- 
inated by the fluctuations in the gravitational potential 
leading to Sachs- Wolfe effect. From FigU we see that 
the Sachs- Wolfe plateau obtained from quasi-exponential 
inflation is almost flat conflrming a nearly scale invari- 
ant spectrum and resonating with the small spectral tilt 
1 — ng — 0.0254 as estimated from our analysis. 

For larger I, CMB anisotropy spectrum is dominated 
by the acoustic oscillations of the baryon-photon fluid 
giving rise to several peaks and troughs in the spec- 
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FIG. 5: Variation of CMB tensor Cf^ angular power spec- 
trum for quasi-exponential inflation and the predicted tensor 
spectra of WMAP7 for ACDM -|- TENS with the multipoles 
I 



FIG. 6: Variation of CMB angular power spectrum Cf for 
quasi-exponential inflation and the best fit spectra of WMAP7 
for ACDM -I- TENS with the multipoles I 
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trum. The heights of the peaks are very susceptible to 
the baryon fraction. Also the peak positions are sensitive 
to the curvature of the space and on the rate of cosmolog- 
ical expansion hence on the dark energy and other forms 
of the matter. In FigU] the first and most prominent 
peak arises at I — 220 at a height of 5811//K^ followed 
by two equal height peaks at I = 537 and I — 815. This 
is in excellent agreement with WMAP seven years data 

for ACDM background. The direct comparison of our 
prediction for p = 0.60 in Fig|4] shows fine match with 
WMAP data apart from the two outliers at Z = 22 and 

1 = 40. The above analysis is for a representative value 
of the parameter p — 0.6 witin its allowed range. 

The gravitational waves generated during inflation also 
remain constant on super Hubble scales having small am- 
plitudes. But as its wavelength becomes smaller than 
horizon the amplitude begins to die off very rapidly. 
So the small scale modes have no impact in the CMB 
anisotropy spectrum only the large scale modes have lit- 
tle contribution and this is obvious from FiglSJ 

Further, in FiglHand FigHwe have plotted CMB TE 
and EE angular power spectrum respectively for quasi- 
exponential inflation and the best fit spectra of WMAP7 
for ACDM -I- TENS and compared with WMAP seven 
years data. Both the plots resonate fairly well with the 
latest WMAP data 

Thus, from the entire analysis, it turns out that 
quasi-exponential inflation confronts extremely well with 
WMAP seven dataset. We expect that the forthcoming 
data from PLANCK would constrain quasi-exponential 
inflation further. 



V. CONCLUSION 

In this article we have confronted quasi-exponential 
models of inflation with WMAP seven year dataset us- 
ing Hamilton- Jacobi formalism. We have first developed 
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FIG. 7; Variation of CMB angular power spectrum Cf*^ for 
quasi-exponential inflation and the best fit spectra of WMAP7 
for ACDM + TENS with the multipoles I 



the formalism with a phenomenological Hubble parame- 
ter and have demonstrated how and to what extent the 
scenario measures deviation from standard exponential 
inflation [e.g., de-Sitter model). The deviation, incor- 
porated through a new parameter p, has then been con- 
strained by estimating the major observable parameters 
from the model and confronting them with WMAP seven 
year dataset. The observable parameters, as obtained 
from our analysis, turn out to be in good fit with the 
latest WMAP data. 

We have further utilized the publicly available code 
CAMB [13 to compare with WMAP seven data Q the 
matter power spectrum as well as CMB angular power 
spectra for the TT, TE and EE modes obtained from 
our analysis. This allows us to put stringent constrains 
on the model parameters. CAMB outputs for a emblem- 
atic value of the parameter p within its allowed range 
0.586 <p< V2 are in excellent agreement with the latest 
WMAP data. Values of the most significant cosmolog- 
ical parameters have also been calculated using CAMB 
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and have found to fair well with the observational bounds 
as given by WMAP seven years data. This leads us to 
conclude that quasi-exponential inflation confronts ex- 
tremely well with WMAP seven within a certain param- 
eter space. 

Nevertheless, another appealing aspect of our analysis 
is the possibility of verifying quasi-exponential inflation 
by PLANCK. The current observational bound on the 
ratio of tensor to scalar amplitudes as given by WMAP 
is r < 0.36 at 95% C.L. With the ongoing satel- 

lite experiment PLANCK promising to detect it down to 
the order of 10~^ (l6j in near future, this parameter is 
now playing a pivotal role (along with non-gaussianity) 
to discriminate among different models. The numerical 
estimation reveals that the parameter r is of the order 
of 10^'^ for quasi-exponential inflation. Thus, the forth- 
coming data from PLANCK can be confronted directly 
with quasi-exponential inflation. We are looking forward 
to confront this type of inflation with PLANCK and to 



have more exciting results in near future. We also plan 
to engage ourselves in estimating cosmological parame- 
ters using another highly useful code COSMOMC [3(| . 
Results in this direction will be reported shortly. 
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